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A SEMI-ANALYTIC APPROACH TO THE 


SELF-INDUCED MOTION OF VORTEX SHEETS 
SUMMARY 

The rolling-up of the trailing vortex sheet produced by a wing 
of finite span is calculated as a series expansion in time. For a 
vorticity distribution corresponding to a wing with cusped tips, 
the shape of the sheet is found by summing the series using Pade 
approximants . The sheet remains analytic for some time but ultimately 
develops an exponential spiral at the tips. The centroid of vorticity 
is conserved to high accuracy. 


IV 



1 . 


Introduction 


Within the context of the potential theory of fluid motion sur- 
faces of velocity discontinuity may be characterized as vortex sheets. 
Typically a shear layer m a real fluid is idealized by collapsing 
the region of large velocity gradient onto a sheet across which the 
magnitude and/or direction of the fluid velocity experiences a finite 
jump. The effects of viscosity are not considered; thus this concen- 
trated vorticity cannot diffuse and the sheet will remain of zero 
thickness for all time. The sheet can, however, deform and stretch 
under the influence of its own induced velocity field. For example 
the vortex tube surrounding the circular jet has been observed to 
"roll up" into periodic spirals. Similar behavior has also been 
observed in the two-dimensional analogue. Rosenhead (1931) intro- 
duced a discrete vortex approximation to study the time evolution of 
a sinusoidally perturbed two-dimensional vortex sheet across which 
the velocity jumps discontmuously . The continuous sheet is replaced 
by a line of point vortices whose strengths vary sinusoidally. The 
induced velocity at a given point vortex, at any instant m time, is 
given by the vector sum of the contributions from all the others. In 
the initial stages of motion, at least, the rolling-up phenomenon is 
clearly indicated m his results. Not all vortex sheets must deform, 
however. Two-dimensional gravity waves may propagate without change 
of form on the interface between two fluids of different constant 
density. Within each fluid the motion is irrotational with a discon- 
tinuity in speed across the interface. The search for a wave of con- 
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stant form may be thought of as the determination of that particular 
vortex sheet configuration which preserves its ’'initial" shape. 

A particular vortex sheet whose evolution has been a topic of 
extensive study is that left in the wake of a lifting surface or wing. 
The span-wise distribution of vorticity m the sheet is produced by 
the shedding of vortex lines as the lift varies along the wing from 
centerline to tip. A classical two-dimensional problem is obtained 
by considering the self-induced motion m a plane so far behind the 
wing that the bound vorticity produces negligible effect (the 
"Tref ftz-plane" ) . Thus Westwater (1936) computed the two-dimensional 
motion of a sheet of finite span, that is initially flat, with the 
variation of vorticity produced by a wing on which the lift varies 
elliptically . Following Rosenhead, he used the discrete vortex 
approximation using 10 point vortices of the same strength distri- 
buted along a semi-span. Each vortex is placed initially at the cen- 
troid of vorticity of the segment of the distribution it replaces. 

His results suggest an orderly rolling-up pattern starting at the edges. 
Earlier, Kaden (1931) found an analytic expression for the form of 
the sheet in the neighborhood of the edge. He considered a semi- 
infinite sheet of vorticity produced by the parabolic lift variation 
that approximates the elliptic variation at the edge. Because there 
is no charactensitc length in the field, Kaden was able to extract 
a simple similarity solution for the shape that, m polar coordinates, 
is a spiral whose form is given approximately by 
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where t is time. Thxs is a tightly-wound spiral of infinite length, 

2/3 

a typical dimension of which grows as t . Kaden's results confirm 
a still earlier prediction of Prandtl (1927) that vortex sheets behind 
wings will roll up towards their tips. Because Kaden's similarity 
solution suggests that the spiral is always completely wound up, even 
at t = 0 + , the sheet, while sensibly flat at t = 0, is, m fact, 
higher singular at the tip. 

With the advent of automatic computing, it became possible to 
pursue the discrete-vortex or "multi-vortex" approximation m much 
greater detail. Takami (1964) and others were unable to reproduce 
Westwater's smooth roll-up results. Chaotic motion was observed near 
the tips even at the early stages of motion. Presumably the smooth 
patterns obtained by Westwater were due to fortuitous inaccuracy m 
his time-integration scheme. Even more disconcerting was the effect 
of increasing the number of discrete vortices. Rather than improve 
the results, the choatic motion was amplified. When two discrete 
vortices moved "too close" they induced inacceptibly large velocities 
upon each other. Takami also considered other distributions of 
vorticity including the one produced by a wing with cusped tips (3/2- 
power loading). In this case, the strength of the sheet goes to zero 
at the tips and smoother behavior might be expected. On the contrary, 
his results for this case indicate that the region of disorderly 
motion is not confined to the vicinity of the tips but extends over 
much of the sheet. 

A serious criticism of the multi-vortex approximation was made 
by Birkhoff & Fisher (1959) . They assert that the self-induced 
motion of an array of point vortices will ultimately produce randomness 
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of position and hence that no true rolling-up is possible. Perhaps 
motivated by this objection, several authors have modified the numer- 
ical problem through the introduction of various smoothing techniques. 
Thus Chonn & Bernard (1973), for example, introduce an arbitrary 
maximum on the permissible induced velocity and claim that this pro- 
cedure reproduces some of the features of viscosity. 

A significant step foward was made by Fink and Soh (1974, 1978) . 
After carefully comparing the multi-vortex model with the original 
Cauchy principal-value integral, they concluded that the former in- 
volves the neglect of initially-small logarithmic terms. This error 
becomes amplified as the sheet moves and ultimately leads to the 
observed chaotic motions. Their improvement is simply to rediscretize 
at each time step. In a number of applications, including the rolling- 
up of a trailing vortex sheet, their results remain smooth for much 
longer periods of time than had been previously reported. 

In addition to its inherent mathematical interest, the problem 
of trailing vortex-sheet roll-up is of significant practical importance. 
Spreiter & Sacks (1951) show that for heavily- loaded low-aspect- 
ratio wings, the sheet may become essentially rolled up into two 
vortex cores within a chord length of the wing trailing edge. This 
effect must be considered in a valid analysis of plane tail perfor- 
mance in these cases. Additional interest followed the introduction 
of wide-bodied transport planes with heavily loaded wings of high 
aspect ratio m the early 1970's. Strong vortex cores left by the 
passage of these large aircraft in the vicinity of airports have 
been implicated as a contributing factor to accidents involving 
smaller aircraft. 
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In the present work we seek to solve for the self-induced motion 
of a finite two-dimensional vortex sheet without introducing any dis- 
cretization at all. The position of the vortex sheet is calculated 
as a power series in time with coefficients that are analytic func- 
tions of a curve parameter. Thus the procedure is restricted to 
those lift distributions that produce sheets whose motion is analytic 
in time initially. 

In §2 we show that for a particular class of vorticity distri- 
butions 1C (x^) , the coefficients m the series will be polynomial 
functions of x^ the curve parameter. For this class an algorithm 
to find the coefficients, suitable for machine computation, is worked 


out in detail. We pay particular attention to the case of 3/2-power 

........ 2, 3 /2 . . _ 


lift distribution, i.e. (1 - x 


m dimensionless units. The 


corresponding vorticity K (x ) = 3x vl - x 2 is the most singular 

o o o 

distribution in the class for which the Holder condition [see e.g. 


Muskhelishvili (1958) ] 


I K (x )| < A I x - l| a , 0 < A < °°, 0<a<l 

o o 

is satisfied at the edge(s) of the sheet. The Holder condition, which 
is stronger than contmuty but weaker than differentiability, is suf- 
ficient to ensure existence of the Cauchy principal-value integrals 
from which the series coefficients are calculated. 

Results of the algorithm developed in §2 are presented m *3 3. 

42 

The series for the sheet coordinates is computed to 0(t ) for the 

3/2-power lift case. The series have a finite radius of convergence; 
convergence fails first at the tips for a dimensionless value of time 
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of about 0.39. Using Pade approximants highly accurate sheet profiles 
are computed and are compared with the numerical work of Takami and 
Fink & Soh. The limiting singularity is associated with the instan- 
taneous appearance of a loosely-wound spiral of finite length. Unlike 
Kaden*s power-low spiral, the local solution here is 



Once the spiral appears at the tip the analytic solution can no 
longer be used. One may conjecture, however, that, as time proceeds, 
vorticity is drawn into the vortex core until, ultimately, it is all 
concentrated there. 

A useful check on the series solution is provided by the invar- 
iance of the horizontal coordinate of the vorticity centroid. Two 
separate checks can be formulated: one relating to a weighted sum 

of the series coefficients at any order and, secondly, a global check 
involving numerical integration of the Pade-summed series results. 

Both suggest that the present results are effectively exact until the 
critical time is reached. 

Finally we compute the series solution that is associated with 
a slightly-perturbed elliptic lift variation. This ’'solution” is com- 
pletely analytic and does not include the singularity at the tip at 
t = 0 . Thus it is incomplete, does not conserve the position of the 
centroid, but it may be useful as an "outer” solution for purposes of 
matching. 
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2 . 


Mathematical Formulation 


The velocity field induced at time t, by a vortex sheet with 
concentrated vorticity distribution K(s,t) where s is arc length, 
is given by 


q(z,t) 


1 f K (s,t)ds 
2 t t J z - z (s , t) 


( 1 ) 


c 

m the usual complex notation. The induced velocity at points on 
the sheet is also given by this expression if the Cauchy principal 

'b 

value of the integral is taken. K is, in fact, equal to the differ- 
ence m the tangential components of velocity across the sheet. 
Assuming that, at the initial instant of time, the sheet lies on 
the x-axis between -b and b we introduce the line parameter x e[-b,b] 


and the "Lagrangian" vorticity distribution k(x q ) defined by 


K (s , t) 


K(X ) 

o 


dx 
o 

ds 


Since the fluid is assumed to be inviscid, the time-dependence 
is only found m the sheet-stretching factor dx^/ds. Equation (1) 
becomes, for points on the sheet. 


u(x ,t) - iv(x ,t) 
o o 



K(€)dg 

z(x ,t) - z(£,t) 
o 


( 2 ) 


The lift distribution is taken to be bilaterally symmetric; hence 

K (x ) is antisymmetric and we introduce the circulation about half 
o 

the vortex sheet 


r = 



K(x )dx 
o o 


( 3 ) 
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The problem may be made dimensionless be selecting b as reference 
2 _ 

length and b /T as reference time. Separately equating real and 

g 

imaginary parts of (2) and recognizing that u(x ,t) = -r— x(x ,t) and 

o d t o 

g 

v (x , t) = 7 rry(x , t), we obtain the coupled system of nonlinear mtegro- 
o d t o 

differential equations 



t) 



ic(£) [y(x Q ,t) - y(5,t)] d£ 


(4a) 


and 


where 



,t) 



K (£) [X(X ,t) - x(£,t) ] 

O 


2 

R 




(4b) 


R = [x (x ,t) 
o 


- x(£ f t)] + 


[y (x ,t) 
O 


- y (C/t) ] 


(4c) 


We shall now seek solutions of (4) which may be developed as power 
series m time that are uniformly convergent for x^t-1,1]. Clearly 
if such solutions exist they must be of the form 


x( V t} = X o + ^ A x (X o )t21 
1=1 


(5a) 


y (x Q ,t) 


-2 


2i+l 

B (x ) t 
l o 


(5b) 
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If we further assume that the vorticity is of the form 


K(S) 



(6) 


it can be shown that the coefficients and will be polynomial 
functions of their arguments. The dimensionless form of equation (3) 
requires that 


c + 
o 


I 


c 

3 


( 2 « 4 » 6 * * * 2 ^ 

3* 5 • 7* • • (2d -t- 1) 


= 1 


The form (6) includes the elliptic lift distribution and the 
3/2-power distribution produced by a cusped wing planform as special 
cases. It does not include the distributions characteristic of rectan- 
gular or lenticular planforms; indeed, no series development in time, 
starting from an initially flat sheet, is possible for these later 
cases. 

When (5) and (6) are substituted in (4) and the coefficients of 
the various series m t are collected, the Cauchy integrations can 
be performed using the family of "airfoil integrals" 


& 

n 



g n dg 

A - 5 2 (x - o 

O 


n — 1 , 2 , 


(7) 


The expressions for & (x ) may be found recursively according to the 

no 

scheme 
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(n even) 


(8) 


& x = -IT 


& - x & - 

n o **n-l 


^ 1* 3* 5* • • (n-2) , 

= X & _ - 7T ■ - ; — — (n odd) 

°n o n-1 2*4*6* • • (n-1) 


It is important to note that are defined only on the open 

interval (-1, 1). 

The series expansion procedure will now be described in detail. 
Several intermediate variables will be introduced both for computa- 
tional convenience and also to reduce the problem to one that is 
only quadratically nonlinear. We define A q = x q and let 

C = A (x ) - A (?), (9) 

1 1 O 1 

and 

D. = B. (x ) - B (£) . (10) 

1 1 O 1 


From equation (4c) , 


2 V ^ 2k 

E - Z V 

k=0 


where 


E = C c +Vd_d 
k i k-i l-l k-i 


1=0 


1=1 


Let 


y(x Q ,t - y(£,t) “ 

~ ^ = ^ n 

3=0 


P t 
1 


2]+l 


( 11 ) 


( 12 ) 


( 13 ) 
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and 


x(x ,t) - x(£,t) 

o 


R 


2 


oo 



D=0 


(14) 


After equating coefficients of like powers of t, we obtain, from 
(13) and (14) respectively. 


E F 
o i 


D -Y 


E k F i~k 

k=l 


(15) 


E G 
o 1 


“ C x-I 


E G 
k l-k 


k=l 


(16) 


for all positive integers 1 . The summations are identically zero 
when i=0. Finally we differentiate equations (5) with respect to 
t and substitute into (4a) and (4b) to obtain 

- 41 T (1 + 1) A 1+1 = j K{£ > )F ± d£ , (17a) 

-1 

2TT(2i + 1) = ^ 1 K(?)G 1 d£ , (17b) 

-1 


for l > 0 . 

Equations 9, 10, 12, and 15-17 form a complete system for the 
successive determination of the coefficients If K (£j) is of 

the form (6) , A. and B will be polynomials m x and C , D , E , 

i i Q ill 

(x - £)F and (x - C)G are polynomials in x and £ . The degree 
o i o i o 
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of these polynomials, for each i, will increase, in general, with 
increasing N. If, for example, N = 1 in equation ( 6 ) it can be 
established by induction that a sufficient form is given by 


W 


-I 


a .x 
1 ] o 


23+1 


D=0 


(18a) 


l+l 

W - I 6 

j=0 


x 

13 o 


and 


where 


2 k 

' (x o ■ 5)2 2 ' 

t =0 


2i+l 


F = — 
i (x 


■W 2 


F (x )£' 
is o 


s =0 


2i 


5 1 " (x/- D 2 


g (x >r , 

is o 


s =0 


kt 


- 1 


E x 2p 
ktp o 


ft +1 

3= [—J 


F 

is 


l+l 


j=l+ [s/ 2 ] 


_ 23-l-s 

F . x J 
IS3 o 


(18b) 


(18c) 


(18d) 


(18e) 


(18f ) 


(18g) 
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and 


G 

is 



G 

is 


2 j-s 

. x 
D o 


(18h) 


Here J denotes the integer-part function. Substituting (18a) , (18b) 

and (18f) into (12) and using identities of the form 

2 d-i 

2j V 2j-l-r >-r 

X / - K = (X - 5 ) > X J 5 

O O O 

r=0 

to eliminate the explicit dependence on the intermediate variables 
and we obtain, after some manipulation. 


E 


ktp 


i 


Min [p,i] 

I 


j=Max [0,p+i-k] 


a a. . y. . 

13 k-i,p -3 tp 3 


♦ 2 


Min Ip, l] 


2 

j=Max [1, p+i-k-1] 


3 3 Y 

i-l ,3 k-i,p- 3 +l T tpj 


(19) 


The new functions introduced m (19) are defined by 
Y t = 1 + Mm [t, 2p-t , 2p-2j ] 
and y = Mm [t+1, 2p-t+l , 2} , 2p-2 j+2] . 

up j 

Summations in (19) are taken to be identically zero when the lower 
limit exceeds the upper. 

Expressions for the triply-subscripted elements an< ^ ^is} 

may be obtained from (15) and (16) as 


13 





2 ]-l-s 

x 


1 Min [2k, s] 


i-k +1 


S S Sr t+ n S r s _ tl E ktp F i-k,s-t,m X i 

k=l t=Max [0,s+2k-2i-l] p= FjM m=l+ [^1 


2 p+ 2 m-l-s 

o 


(s = 0 , 1 , . . . 2 i + 1 ) 


( 20 ) 


and 




lSJ 


a x 2;| - s 
13 / o 


i Min [2k, s] 


i-k 


■- 1 i 


2 2 


^ ^rt+n Ms-t+r 

k=l t=Max[o,s+2k-2i] p= I- m= ■ — 


E G. x 2p+2m s 

ktp i-k,s-t,m o 


(s = 0 , 1 , . . . 2 i) . 


( 21 ) 


Note that the last subscript in F and G on the left- 

133 1S3 

hand side of ( 20 ) and ( 21 ) is given implicitly or computed from the 
exponents of x q in the quadruple sums. While explicit expressions 
similar to those obtained for E^ t ln (19) are possible in principle, 
the present procedure is at least as efficient computationally and 
avoids a great deal of laborious analysis . 

The system of equations for the coefficients is complete once 
the reduced form of equations (17) is obtained. For the special case 
of 3/2-power wing loading, c q = -c^ = 3, the induced velocities are 
continuous at the edges of the sheet and equations (17) can be written 


14 



in a relatively simple form. The relevant family of Cauchy integrals 


is 


i = — <£ 
n t r J 


-1 




dC 


n — 1 / 2, 


which are related to the original family according to 


*n 7T ( ^n ^n+2 , . 


I n can be computed recursively as 


I, - ** - X /2 - 


I = X I , 

n o n-1 


(n even) , 


I = x I 
n o n -l 


1*3*5* * * (n-2) 
2*4*6* • * (n+1) 


(n odd) , 


from equation (8) . 

They are even or odd polynomials m x , of the form 

o 


m 

■■-2 


K x 
p o 


n+l-2p 


p=0 


where K are determined by 
P 


K = 1 
o 


K = 2P ~ L K 
p+1 2p+2 p 


p = 0, 1. 


( 22 ) 
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Equations (17) now become 


l+l 


^ a i+l,£ 

1=0 


21+1 

X 

o 


2i+l 

j&r I 

s=0 


l+l 




K F. 

P xsj 


x 2(n-p)+l 

o 


(23) 


and 


l+l 



3 

2 (2i+l) 


2 

3=n 



(24) 


Starting with = 1 and i = 0, the five coefficient equa- 
tions are solved in the order (19), (21), (24), (20), (23). The index 

l is then incremented by one and the next order of calculation is 
performed. 
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3. 


Discussion of Results 


The special case corresponding to the 3/2-power lift distribu- 

42 

tion has been computed to 0(t ) using an optimized F0RTRAN compiler 

on the Stanford IBM 370. Execution time for a solution of order 

2N 5 

t was found to be proportional to N which is consistent with 

the number of nested summations m the algorithm of the last section. 

For 2N = 42, the computation required 2.2 sec of CPU time. Double- 

precision (16-place) arithmetic was used and all coefficients are 

accurate to at least 4 places even at the highest order. Extended 

precision would be necessary however for a solution of still higher 

order. 

7 

Though 0 (t ) the coefficients m the double series for x and 
y can be recognized as rational numbers: 


x(x ,t) 
o 


x 

o 


+ X 

o 



9 

4 



2 


t 


135 _ 135 2 27_ 4 

128 32 X o 8 X o 


+ x 

o 


8343 

5120 


31509 2 4131 4 

2560 X o 160 X o 


81 6 \ 6 „ ,8, 
— x It + 0 (t ) , 
5 o J 


(25a) 


y(x ,t) 

O 





81 2 
128 X o 


81 4 567 6 \ 5 

— ■ X + T— r X t 

32 o 160 o J 


(25b) 


+ 


7047 

14336 


243 2 

112 X o 


73629 4 7533 6 

17920 X o 320 X o 


90639 8 

4480 X o 


t 


7 


9 

0(t ) 
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4 

The coefficients in (25) through 0(t ) agree with those calculated 

by Professor M. D. Van Dyke as reported by Takami (1964) . 

For a given value of t, the coordinates of a point (x,y) lying 

on the vortex sheet can be found by summing the series m equations 

(5a) and (5b) . Since only a finite number of terms m these series 

are known, their sum, for each value of x q , can be approximated by 

considering the convergence of the sequence of "diagonal Fade 

approximants" formed from such series. Pade approximants are ratios 

of polynomials with coefficients so chosen so that, when expanded 

for small argument, the power series expansions of these ratios 

agree with the original series to appropriate order. Diagonal 

approximants have the additional property that the order of the 

numerator and denominator are equal. Thus if the power series for 

2N 

a function f (G) is known through 0(G ), 

f(e) = a + a. £+...+ a_ g^ n + 0(g^ n+ ^) 
o 1 2N 


we can, in general, form an approximant, denoted by [N/N] f , so that 


f (e) = [n/n] f + 0(e 2N+1 ) , 


where 


[N/N] f = 


b + b G + . . . + b G 

o 1 N 


N 


1 + C,G + . . . + C G 
1 N 


N 


and the b and c are determined uniquely from a , . • - , a_ T for 
3 3 o 2N 

given N. The sequence of approximants so formed will usually con- 
verge much faster than the original series and can converge to the 
analytic continuation of the series if G lies outside the radius 
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of convergence. A number of examples of the use of this technique 
m the solution of fluid mechanics problems can be found m Cabannes 
(1976) . 

Figure 1 shows two configurations of the vortex sheet drawn for 

2 

t - .05 and .15. These have been computed as [10/10] approximants 
2 

to the t -series for x(x ,t) and y(x , t)/t. Also shown in the 

o o 

figure is the induced velocity profile at t=0 m arbitrary units. 

2 

For t = .05, the sequence of diagonal approximants converged to 

14-decimal places for all x^ in [0,1]. Even the [3/3] approximants 

constructed from only seven terms m each series in (5) agree to 

2 

7-decimal places with the converged results. For t = .15, on the 

other hand, 5-place convergence was obtained for x^ <_ 0.998. We 

2 

will show below that a singularity has appeared in the t -series 

to destroy convergence for x slightly greater than this value. 

2 

Notice that the sheet has curved back for t = .15; the value of 
the parameter x q where the tangent to the sheet is vertical is 
about 0.975. 

In Figure 2 we compare the results of the present method with 

the multi-vortex results of Takami (1964) and the numerical results 

of Fink and Soh (1974) . Both authors have produced solutions for 
2 

t = .16. Takami ' s results exhibit the partial disorder character- 
istic of the simple multi-vortex representation. Fink and Soh have 
computed a smooth shape for the sheet that is in general agreement 
with the [10/10] approximant solution except near the tip of the 
sheet. Their relatively coarse point spacing is, apparently, unable 
to resolve the details m this region. In fact, the sheet exhibits 
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2 

a singularity at x^ £ 0.997 for this value of t . Convergence of 
the diagonal approximants failed for x q > .990 and we have only 
drawn the sheet up to this point. 

For a bilaterally symmetric vortex sheet it can be shown that 
the horizontal position of the centroid of vorticity, which m our 
parametric notation is defined by 


X = 



0 


x(x ,t)K(x )dx 
o o o 


is an invariant of the motion. We have 


(26) 


dX 

dt 


= f |jr(x ,t)K(x ) dx 
J dt O O O 

0 


(27) 


Using equation (4a) , the right side of (27) can be written as 



k(x ) dx 
o o 


f 

-1 


1 

k(5) 


[y (x Q , t) - y (£,t) ] 




_l_ 

27T 


jY 

0 0 


[y (x / 1) - y (£, t) ] 

K (x )K(£) d£ dx 

o 2 o 

R 


(28) 



0 0 


K(x )k(£) [ y(x ,t) - y(£,t)] d£ dx 
o o o 

[x (x ,t) + x(£rt)] 2 + [y (x ,t) - y(£,t)] 2 
o o 
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where we have used the symmetry requirements k(-£) = -k(£), y(-£,t) = 
y(£,t), x(-£,t) = -x(£,t) to obtain the last integral. Each integral 
on the right of (28) is invariant under the interchange of the dummy 
arguments x^ and £ and hence is equal to zero. 

The fact that dX/dt is zero can be used in two different ways 
to check our solution. Substituting expansion (5a) in (27) we obtain 
immediately 


/ 


1 


A. (x )K (x ) dx 
1 o o o 


= 0 , 


0 


i 1 / 2 f . . . 


V 2n+l / 2 

With A. = > a x and K = 3x vl-x f the integrations can 

l Z-f 13 o 00 

D“0 

be performed to yield a check on the sums of the coefficients at each 
order in is 


where 



0, 1 = 1, 2 , . . . 


1*3*5** * ( 2 j-f 1) 
j 2 • 4 * 6* • • (2j+4) 


Performing this check numerically on a produces a result for 
each 1 that is at least 16 orders of magnitude less than the largest 
coefficients m the sum. Thus the solution satisfies the consistency 
relation at each order; we also have an estimate for the round-off 
error m the coefficients. 


21 



The invariance of X can also be used to check the Pade-summed 


results for x(x ,t) . Using 160 equally spaced values of x on 
o o 

2 

[0,1] and Simpson's rule to integrate (26) numerically for t = .14 
produced a result which differed by only 0.02 per cent from the correct 
value X = 37T/16. 

For the 3/2-power lift distribution that we have treated, the 

vortex sheet is an analytic curve for [x q | _< 1 m the initial 

stages of motion. A singularity, that is always present on the 

analytic continuation of the sheet, i.e. |x | > 1, moves inward as 

2 2 

time progresses and arrives at x = 1 when t = t* . 

2 

The nature of the singularity and the corresponding value of t* 
can be estimated by use of a graphical procedure due to Domb & 

Sykes (1957). They note that if 

00 

f(e) = ^ a^ 11 = k(£*~ e) a , a / o, 1, . . . (29a) 

n=o 

then 

a/a . = e" 1 [1 - (1 + a)/n] (29b) 

n-1 * 

which follows from the binomial expansion. Thus, for these special 

cases, if we plot the ratios a /a _ versus 1/n , the points will 

n n-1 

be on a straight line. In general the unknown function f can be 
thought of as the sum of a number of singularities; if the one closest 
to the origin of E is of the above type, then the "Domb-Sykes plot" 
will ultimately tend towards a straight line as 1/n becomes small. 
Values of £* and a appropriate to this singularity can be found 
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from the inclination and intercept of this straight line. When only 

a finite number of ratios are know, estimates of e* and a can 

still be made if the plotted points tend towards a straight-line 

asymptote. It has not generally been recognized that (29) is 

still valid for complex values of a n , m which case the asymptote 

will be a straight line in (l/n,,3??e{a /a n } /t ^m {a /a _})-space. 

n n-1 n n-1 

2 

At the tip of the sheet, x =1, the coefficients m the t - 

o 

series for x and y in equations (5) have been used to construct 

the ratios A (1)/A _ (1) and B (1)/B _ (1). They are given in Table I. 

While each set of ratios appears to tend towards a limit, a careful 

examination reveals that these limits, if they exist, are somewhat 

different. Because of the intimate coupling between the two series 

implied by the governing equations (4a, b) , this difference m the 

limits is unacceptible . That is, for a given value of x q , the series 

2 

for x cannot converge for certain values of t while the series 
for y, which is derived from it, diverges for these values. A more 
consistent interpretation is obtained by considering the complex 
series for the quantity 

x(l,t) + i y(l,t)/t = [A (1) + iB (l)]t 2n = C (l)t 2n . (30) 

/ , n n n 

n=o n=o 

The real and imaginary parts of the ratio (1) / c n _^ (D are 

given m Table II and are also plotted in Figure 3. The abscissa has 

been taken as l/(n + 5/4) ; the shift m n, related to the local 

behavior of a regular function which multiplies the singular one near 
2 

t* is determined so as to minimize the curvature m the real-ratio 
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plot. The straight line shown m the figure has the equation 






- 

c /c , 
n n-1 

_ 1 

i i - 5/2 

f 6.7286 

n + 5/4 

l J 


_ — 


This line provides a virtually perfect fit to the ratios for n >_ 7. 

The points corresponding to the small imaginary part of the ratios 

do not lie on a straight line. They can be fitted to a smooth curve, 

however, which plausibly can be extrapolated to the origin. Near 

2 

the critical value t* = 0.14862 + i 0.0, the left side of (30) varies 
locally like 


_ 3/2-. 067i 


R^t 2 ) + R 2 (t 2 ) 




(31) 


The imaginary part of the exponent is estimated from the slope of the 

curve fitting *^m{c /C J at the origin. Consequently, its value 

n n— 1 

connot be considered as accurate as ^Pe(a) = 3/2 . Because of the 

shift m the horizontal axis in Figure 3, the regular function 

-5/2 

probably behaves locally as t . The regular functions R^ and R 2 
could be estimated by comparing the model equation (31) with the 
original series but this has not been attempted. 

From the model function (31) the local behavior of the vortex 
sheet, specifically the trajectory of the sheet tip, may be deduced. 
Let 


Z 


Re 


0-0 ) 
o 



A J^x ( 1 , t ) - 
3/2 r 


exp 


x (1, t ) + 1 


{y U,t) - y (l,t ) } 
* 
t 


. 067i £n 



(32) 
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where A and 0 q are real constants, representing the dilatation 
and rotation necessary to match the local behavior to the "outer" 
solution. Equation (32) represents a spiral of the form 

2 /3 

(0 - 0 ) = - .067i in R 

o 

or 

-22.4(9-0 ) 

R = e ° (33) 

Unlike Kaden's (1931) similarity solution for the tip region of 
the vortex sheet shed by an elliptically-loaded wing, R = K/0 
which is a tightly-round spiral of infinite length, (33) represents 
a loosely-wound spiral whose length, as 0 varies from 0^ to 
infinity, is finite. The numerical constant m (33) shows that in 
one revolution the spiral radius decreases by about 60 orders of 
magnitude! 

Taking the time derivative of (31) reveals that the velocity 

* 

components remain finite as t t at the sheet tip. Since the 

sheet is continuous for all t < t* , the spiral m (33) can also 

represent the spacial form of the sheet m the neighborhood of the 
* 

tip at t = t . Because of the exponential character of the spiral 

and its very loose winding, it is not surprising that the roll-up 

cannot be observed m the sheet configurations m Figures 1 and 2. 

As mentioned above, the radius of convergence of the series m 

2 

time varies with x q . The procedure used to predict t^ for x q = 1 
was also employed for x^ = 0.999 and 1.001 to yield the first-order 
variation 

t 2 (x ) = .14862 - 2.88 (x - 1) 

* o o 
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near x =1. For significantly smaller values of x , the ratio- 
o o 

plot procedure did not give accurate results, presumably due to the 

complexity of the pattern of singularities. It is clear, however, 

that the radius of convergence increases as |x q | is reduced. For 
2 

values of t greater than ,14862, the series solution cannot predict 

the evolution of the vortex sheet because the vorticity distribution 

is no longer analytic on it. One may conjecture however that the 

vorticity between the critical value of x and x = 1 becomes con- 

o o 

centrated at the center of the spiral. The amount of vorticity at this 
point will increase with time; if, ultimately, all the vorticity be- 
comes concentrated there, this point must lie at the centroid location, 
X = 377/16, and will move downward with the constant speed predicted 
for a counter-rotating vortex pair. 

As the sheet deforms from its initially flat configuration, 
vorticity is convected outward towards the tips. The vortex intensity 
is given by 

K(s,t) = K(x )/(ds/dx ) (34) 

o o 

which varies as the sheet is stretched. Since 



the stretching may be computed using Pade sums for the series 
2 V* 2n 

L A ' (x )t and 1 2* B 7 (x )t . Figure 4 shows dx/dx plotted 
no n n o o 

n 2 

versus x for t = 0, 0.05 and 0.148. Combining these results 
o 

'X, 

with K (x ) gives K(s,t) according to (34). It is plotted versus 
o 
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'V 

x in Figure 5. Note that K is a double-valued function of x near 
2 

the tip for t = 0.148 because the sheet has bent back towards the 
center. 

We shall now consider formal power series solutions for vortex 
sheets whose vorticity distributions differ slightly from that pro- 
duced by elliptic loading. For strictly elliptic loading, the air- 
foil integral (7) predicts constant downwash velocity on the open 

interval x e(-l, 1). To the extent that we restrict consideration 
o 

here to analytic sheet configurations, the "formal" series solution 
is simply 


x(x o' t} = X o ' 

y ( v t)=_ 2 t • 


(35) 


This "solution" is incomplete, however. As x q 1 + , for example, 
infinite upwash results. Because the induced velocities are not con- 
tinuous at the tips, the vortex sheet will not be analytic there. 
Because the discontinuity is infinite, moreover, the vertical induced 
velocity at t = 0 must include, at leading order, a singularity of 
the nature of a Dirac 5-f unction there. This infinite velocity, 
directed at right angles to the tangent to the sheet, causes the tip 
to roll up mtantaneously into the similarity form predicted by Kaden. 
The solution (35) is valid however at sufficient distance from the 
tips for sufficiently small time. It is, in essence, an "outer solu- 
tion" which must be joined in some way to Kaden 1 s description of the, 
initially small, vortex core. 
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Similar "outer" solutions can be produced by the present method 
for slightly perturbed elliptic distributions. Using the algorithm 
described in §2, we consider the two distributions 


and 


<i ( Xo ) 


15x / 2 

14/I^7 \ " 10 

o ' 


K,(X ) 
2 O 


15x 


16/l-x£ 



(36a) 


(36b) 


The series method generates the terms forced by these vorticity dis- 
tributions using the integrals (7); the resulting solution is, of 
course, incomplete but fills the same role as does (35) for the pure 
elliptic case. 

Figure 6 shows the shape of the sheet, computed by the Pade 
method, at a later time, for each distribution m (36). In both cases, 
the time series possess finite radii of convergence; again conver- 
gence fails first at the tips. The horizontal location of the cen- 
troid is not invariant for these incomplete solutions. 
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1 


TABLE 


■ * 

Ratios of Coefficients 

from Equations (5) 

n 

A (1)/A ,(D 

n n-± 

B (1) /B (1) 

n n-l 

1 

-1.125 

-1.125 

2 

.1875 

.825 

3 

5.025 

3.5601 

4 

2.6980 

3.4615 

5 

3.6734 

4.0983 

6 

3.9945 

4.4265 

7 

4.3389 

4.7077 

8 

4.5869 

4.9226 

9 

4.7915 

5.0976 

10 

4.9584 

5.2414 

11 

5.0979 

5.3620 

12 

5.2161 

5.4646 

13 

5.3176 

5.5528 

14 

5.4055 

5.6295 

15 

5.4826 

5.6968 

16 

5.5505 

5.7563 

17 

5.6110 

5.8093 

18 

5.6650 

5.8569 

19 

5.7137 

5.8997 

20 

5.7577 

5.9386 
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TABLE II: Ratios of Complex Coefficients C n (l)/C n ^ (1) from Equation (30) 


n 

3fe { C n ,1>/C n-l‘ 1 >} 

t^mic (1)/C 

l n ] 

1 

-1.125 

0.0 

2 

.4170 

.3060 

3 

3.6833 

-.4066 

4 

3.3434 

.2761 

5 

4.0558 

.1275 

6 

4.3911 

.1185 

7 

4.6827 

.0927 

8 

4.9031 

.0786 

9 

5.0820 

.0673 

10 

5.2286 

.0588 

11 

5.3513 

.0522 

12 

5.4554 

.0468 

13 

5.5448 

.0425 

14 

5.6225 

.0388 

15 

5.6906 

.0358 

16 

5.7508 

.0331 

17 

5.8044 

.0308 

18 

5.8524 

.0289 

19 

5.8957 

.0271 

20 

5.9349 

.0256 
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k = ^(Xq), t = 2.8 
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Figure 6. Sheet configurations from "incomplete" solution for 


perturbed elliptic distributions. 
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